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and control performance limitations 
in subcritical plane Poiseuille flow 
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Suppression of the transient energy growth in subcritical plane Poiseuille flow via feedback control 
is addressed. It is assumed that the time derivative of any of the velocity components can be 
imposed at the walls as control input, and that full-state information is available. We show that it is 
impossible to design a linear state-feedback controller that leads to a closed-loop flow system without 
transient energy growth. In a subsequent step, full-state feedback controllers - directly targeting the 
transient growth mechanism - are designed, using a procedure based on a Linear Matrix Inequalities 
approach. The performance of such controllers is analyzed first in the linear case, where comparison 
to previously proposed linear-quadratic optimal controllers is made; further, transition thresholds 
are evaluated via Direct Numerical Simulations of the controlled three-dimensional Poiseuille flow 
against different initial conditions of physical interest, employing different velocity components as 
wall actuation. The present controllers are effective in increasing the transition thresholds in closed 
loop, with varying degree of performance depending on the initial condition and the actuation 
component employed. 



I. INTRODUCTION 

Transient energy growth has been recognized as a pos- 
sible mechanism explaining subcritical transition in wall- 
bounded flows; in fact, subcritical flows may experience 
large transient amplifications of the energy of perturba- 
tions, that could trigger nonlinear mechanisms and even- 
tually lead to transition to turbulence^"— . 

In viscous shear flows, transient energy growth is re- 
lated to the non-normality of the linearized Navier-Stokes 
operator with respect to the energy inner producf^'^. In 
the last few years, several investigators attempted to re- 
duce the transient growth phenomenon in Poiseuille and 
boundary layer flows by employing wall actuation and ap- 
plying linear control theory to an appropriate discretiza- 
tion of the linearized equations. In their seminal work 
on feedback control of instabilities in two-dimensional 
Poiseuille flow, Joshi et al,® employed a compensator in 
the form of a constant-gain integral feedback, and demon- 
strated stabilization of the linearly unstable flow as well 
as attenuation of finite amplitude disturbances. They 
further pointed out that transient amplifications in the 
flow energy may not be properly detected by the sen- 
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sors, and also that the control itself may trigger non- 
linear mechanisms by introducing transient disturbances 
on short times. Leveraging a state-space formulation 
obtained after discretization of the boundary-controlled 
Orr-Sommerfeld-Squire equations, optimal and robust 
control theory was applied to transitional channel flows 
by Bewley & Liu^ for a single wavenumber pair and by 
Hogberg et al.* for a large array of wavenumber pairs, 
leading to a reduction of the maximum transient energy 
growth as well as an increase in transition thresholds. It 
has been recently shown^ that the linear coupling term 
in the Orr-Sommerfcld-Squire equations plays a role not 
only in the non-normal behavior of the small perturba- 
tion dynamics, but also in the self-sustaining process of 
near-wall, low Reynolds number turbulence. This evi- 
dence led investigators to test in turbulent channel flows 
the optimal controllers designed on linearized flow mod- 
els, and encouraging results have been obtained in terms 
of drag reductioniSrJ^. Feedback control of non-modal 
disturbances in boundary layer flows has been recently 
considered by Corbett & Bottaroi^ in the framework of 
optimal control theory, while Zuccher et al.~ applied 
steady suction in the attenuation of the growth of given 
optimal disturbances in a Blasius boundary layer. 

Although it has been demonstrated^"^ that optimal and 
robust control laws are well suited for reducing the non- 
normal behavior of fluid flow systems, to date no feed- 
back control law has been devised with the capability of 
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ensuring closed-loop monotonic stability, when bound- 
ary actuation is employed. It is therefore natural to ask 
whether such performance can be obtained with a linear 
feedback, and further which control techniques are avail- 
able to directly target the transient growth mechanism. 
In the controls literature, the transient amplification of 
certain norms of the state bears the name of peaking 
phenomenon and the monotonic stability requirement is 
generally referred to as strict dissipativity. Active con- 
trollers with the capability of targeting transient norm 
amplifications have received attention in the analysis of 
a class of partially linear cascade systemsi^, and more 
recently in conjunction with a Linear Matrix Inequality 
(LMI) approacbi^^. In a very recent paper, Whidborne 
& McKernani^ extended these results giving conditions 
on the existence of a feedback controller ensuring the 
strict dissipativity of the closed-loop system. These re- 
sults were then exploited by Whidborne et. al^, who 
considered the feedback control of a single wavenumber 
pair in plane Poiseuille flow via LMI design and wall- 
normal blowing and suction, and by Martinelli et alj^, 
where LMI-based feedback controllers have been designed 
and tested for an array of wavenumber pairs. The present 
paper aims at expanding and completing these recent re- 
sults, by showing first that it is impossible to design a lin- 
ear state-feedback controller ensuring the plane Poiseuille 
flow - controlled via wall transpiration with any veloc- 
ity component - to be strictly dissipative. In a second 
step, feedback control laws are designed using an LMI 
technique, for an array of wavenumber pairs; their per- 
formance is compared against that of optimal controllers 
in the linear case, and furthermore closed-loop transition 
thresholds are evaluated for optimal initial conditions in 
the form of a pair of oblique waves and antisymmetric 
streamwise vortices, at the Reynolds number Re = 2000, 
using different velocity components as wall actuators. 



II. MODEL OF THE SYSTEM 

We consider the dynamics of three-dimensional small 
perturbations to the laminar Poiseuille solution in a plane 
channel. A Cartesian coordinate system is introduced, 
where x, y and z denote the streamwise, wall-normal 
and spanwise directions, and u, v, w denote the corre- 
sponding perturbation velocity components. The Navier- 
Stokes equations, linearized about the laminar solution 
Uiy) = Up{l — (y/cJ)^), are non-dimensionalized with 
the center line velocity Up and the channel half- width 5, 
and rewritten in the form of a single equation for v one- 
way coupled to an equation for the wall-normal vorticity 
rj — du/dz — dw/dx. Fourier transformation in x and 
z direction yields the well known Orr-Sommerfeld-Squire 
form: 

Aw = [-jaC/A + jaU" + AA/i?e]w 



at the wavenumber pair (a,/3). Here, the tilde denotes 
Fourier coefficients, the dot denotes time derivative, the 
prime denotes y differentiation, = + 0^ , j is \/— T, 
and A = (f/dy"^ - k^. 

We select boundary conditions representing time- 
varying wall transpiration on any of the velocity com- 
ponents at the two channel walls ( "vectorized transpira- 
tion"). In turn, this results in inhomogeneous Dirichlet 
and Neumann conditions on v, as well as inhomogeneous 
Dirichlet conditions on fj: 

v{y ±l,t) = Vuj{t), 
dv/dyiy^±l,t)^iUt), (2) 

viy = ±1,0 = Vu,i{t)- 

A standard lifting procedureS is then employed, i.e. the 
unknowns are rewritten as homogeneous components 
(satisfying homogeneous boundary conditions) plus in- 
homogeneous components thus: 

v{y,t) = ih{y,t) + fu{v)iu{t) + Ii{v)ii{t) + . . . 

■ ■ ■ + gu{y)v^{t) + gi{v)vi{t) (3) 

fl{y,t) = f]h{y,t) + K{y)f]u{t) + hi{y)f]i{t), 

where fu, fi, gu, gi, hu, hi are polynomials in y chosen to 
satisfy unitary boundary conditions for each lifted com- 
ponent appropriately. 

We discretise the homogeneous components Vh and 77^ 
in wall-normal direction using a modified Chebyshev se- 
ries cardinal function basis 

N-4 

v,{y,t)^J2^n''iyMt) 

n=0 

Ar-2 ^ ' 

Vh{y,t) = r^{y)ar,,n{t) 

where the modified Chebyshev functions T^^{y) and 
(y) implicitly enforce the required homogeneous 
Dirichlet and Neumann boundary conditions, lead to 
good conditioning of the discrete Laplacian operator, 
and no spurious modes are generated^^. The Orr- 
Sommerfeld-Squire equations are then evaluated on a set 
of Gauss-Lobatto collocation points in y direction and 
rearranged to have the time rate of change of actuation 
velocity as an inpuli^. This results in the the linear time- 
invariant system model^^ 

±(t) = Ax(t) -I- Bu{t), x(0) = xo, 

where A, B are constant system and input matrices, and 
u, X are respective input and state vectors. As noted in 
previous work^, the natural outcome of this procedure 
is an augmented state-space form, where additional in- 
tegrators associated to the values of the input velocity 
components at the walls are explicitly introduced. Con- 
sequently, the open-loop dynamics of the system above 
(i.e. setting _B = 0) is different from the dynamics of the 
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original, unactuated system (where the velocity at the 
walls is fixed by the no-slip condition). This terminology 
will be used throughout the paper to distinguish between 
the two cases. 

The kinetic energy per unit mass of flow perturbations 
in the volume V 

E = ^ I u^+v^ +w^dV 

Jy 

can be expressed as a function of the state vector x using 
the continuity equation, the definition of rj, and Parse- 
val's identity: 

where the matrix Q is an Hermitian, positive definite 
matrix; here, the superscript H denotes conjugate trans- 
pose. In the following, we have transformed the state 
vector via the change of variable x = Cx (C being the 
Cholesky factor of Q) in order to rewrite the system dy- 
namics as 

±{t) = Ax(t) + Bu{t), x(0) = xo, (5) 

such that the system energy is directly given by the Eu- 
clidean norm x^x^i^^. 



III. CLOSED-LOOP MONOTONIC STABILITY 

We consider the linear time-invariant system model 
(HJ and further assume that B^B > 0, that is B has 
full column rank (i.e. all the actuators are independent 
- a condition which is trivially satisfied in the present 
problem). Contraction analysis for this kind of system 
has been presented by Whidborne & McKernaniS and, 
in the general case of nonlinear systems, by Lohmiller 
& Slotine^"*. In particular, referring to the linear case, it 
has been shown^^ that there exists a static state-feedback 
controller u = Kx, where K is a constant matrix, such 
that the closed-loop system has strict dissipativity (i.e. 
energy x^x decays monotonically from all initial condi- 
tions Xo), if and only if 

B-^ (^A + A^) B-^^ < or BB^ > 0, (6) 

where B"*^ is the left null space of B. Additionally^^ if no 
static controller that achieves strict dissipativity exists, 
then no dynamic state-feedback controller, where u is 
given from x by the dynamic system 

ife(t) = AfeXfc(i) + Bfex(t), Xfe(0)=Xfeo, (7) 
u(t) = CfcXfe(t) + Dfex(t), (8) 

where Afe,Bfc,Cfe and are constant matrices and x^ 
are controller states, exists either. 

It is immediate to verify that the second criterion in ^ 
is never satisfied in the present system, as the Hermitian 



matrix BB is never positive definite but it is always 
positive semidefinite, because the dimension of the input 
vector is always smaller than that of the state vector. In 
order to have BB^ > 0, a number of independent actu- 
ators equal to the number of flow states is required; this 
is a situation that is unlikely to occur in practical flow 
control problems, where normally actuators are placed 
at the walls. Even when volume forcing is available, this 
condition is unlikely to be satisfied, since practical vol- 
ume forces are not as flexible as to enforce an arbitrary 
force distribution in the entire flow domain at any time 
instant. 

The first algebraic criterion in Eq. © is equivalent 
to requiring that the portion of the system dynamics 
which is not accessible by the controls must be dissi- 
pative. Verifying this criterion is not trivial. Here we 
evaluate it numerically, in order to identify those re- 
gions in the (a, /?, Re) parametric space where subcrit- 
ical Poiseuille flow may be rendered monotonically sta- 
ble by feedback transpiration. To this aim, the state- 
space model ^ is flrst obtained on a fine grid and, as 
suggested by Reddy & Henningson^, a limited number 
Nt of eigenfunctions is retained, discarding those corre- 
sponding to highly damped and poorly resolved eigenval- 
ues. Properly rescaling the variables such that energy is 
written as an Euclidean norm leads to a reduced order 
model Ar,Br, and the negative-definiteness of the corre- 
sponding matrix B^ ^A^ -I- A^^ B^^ in ([6]) is verified 

by computing its maximum (real) eigenvalue Xmax- Fig- 
ure [T] shows the present result on the (i?e, a) plane for 
/3 = 0, along with the well-known result on the transient 
growth dependence in plane Poiseuille flow^ (i.e. the un- 
actuated case). The white area corresponds to the do- 
main where the unactuated system is monotonically sta- 
ble, while the shaded area is the region where the unactu- 
ated system admits transient energy growth. Solid lines 
correspond to isocontours of Xmax, and it appears that 
the contour Xmax — lies on the very boundary between 
the shaded and white area, implying that the Hermitian 

matrix B^ (^Ar + A^^ B^^ is indeflnite when the unac- 
tuated system is not monotonically stable. An analogous 
result is reported in Fig. [21 where isocontours of Xmax 
are reported on the plane {a, (3), at Re = 120; again, 
the contour Xmax — lies on the boundary between the 
regions of monotonic and non-monotonic stability. From 
the aforestated theorem^^, this implies that it is not pos- 
sible to design a state-feedback controller that ensures 
the closed loop Poiseuille flow to be monotonically sta- 
ble, when the corresponding unactuated flow is not. 

This result shows an inherent limitation in the feed- 
back control of the transient growth mechanism, when 
vectorized wall transpiration in terms of time rate of 
change of zero net-mass flux blowing/suction is em- 
ployed. Note that vectorized transpiration, although be- 
ing rather idealized, exploits all the degrees of freedom 
available for boundary control in the present problem; 
therefore, the present result is representative of a limiting 
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situation. In a practical setting, actuator-dependent con- 
straints may introduce additional mechanisms restricting 
the control authority even more. It is also worth men- 
tioning that strict dissipativity is a rather tough require- 
ment for a controller bound to operate on a largely un- 
deractuated flow, and in fact an approach based on strict 
dissipativity is in general quite conservative, i.e. some en- 
ergy growth is tolerable in transition control. Finally, it 
should be emphasized that the present analysis is limited 
to linear feedback laws, and that performance of nonlin- 
ear controllers may be more promising. For example, it 
has been shown^^ that introducing a nonlinearity in the 
form of gain-scheduling on full-state feedback laws led 
to relaminarization of low-i?e turbulence even employing 
wall actuation only, and that adjoint-based optimization 
on the nonlinear turbulent flow can be successfully em- 
ployed in feedback relaminarizatioi*ii. 
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FIG. 1: Numerical verification of the first algebraic crite- 
rion in Eq. ((6|. Lines: contours at constant Xmax{Re,a), 
at /3 = 0. Levels are from —0.1 to 0.4 by 0.1 increments; 
dashed line is negative value. The shaded area corresponds 
to the region where the unactuated system is not monotoni- 
cally stable, i.e. admits transient energy growth. The contour 
\max{Re,a) = lies on the boundary of the region, indicat- 
ing that no state-feedback controller can be designed to ensure 
strict dissipativity of the closed-loop system when the unactu- 
ated system is not strictly dissipative. Results obtained with 
N = 100, Nt = 50. 



IV. UPPER-BOUND MINIMIZING FEEDBACK 
CONTROLLER 

In order to design a state-feedback controller with the 
capability of targeting the transient growth mechanism 
directly, an estimate of the maximum transient growth is 
required. Such estimate is obtained as an upper bound 
on the maximum growth via Lyapunov theory. For the 
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FIG. 2: Numerical verification of the first algebraic criterion 
in Eq. Results for Re = 120, levels are from -0.1 to 0.2 
by 0.05 increments. For details see caption of Fig. [1] 

linear, time invariant, asymptotically stable system: 

X = Ax, x(0) = xo, 

it can be shown that an upper bound on the maximum 
transient growth G is given b y^^i^° : 

Gu = Amaa; (P)Ama2; (P ^) ^ G, 

where P = P^ > satysfles the Lyapunov inequality 

PA -I- A"P < 0. 

A minimal upper bound can be obtained by solving the 
following minimization problem^^i : 

min7 : 

PA + A^P < 0, P = P^ > (9) 
I< P < 7I, 

where the last inequality ensures 7 > Gu- The problem 
stated in Eq. © is a LMI generalized eigenvalue prob- 
lem, and standard solution methods based on interior 
point algorithms are available^^. 

An analogous problem to that stated in Eq. ^ can 
be obtained if the feedback minimization of the upper 
bound is of interest. Indeed, let us consider the system 
([5]) along with a state-feedback control law in the form 
u = Kx, so that in closed-loop the system dynamics is 
described by: 

x=(A + BK)x, x(0)==xo. 

Leveraging the additional degrees of freedom due to the 
controller gains K, we move to minimizing the closed- 
loop upper bound on the maximum transient growth. 
The associated Lyapunov inequality now reads: 

PA + A^P + PBK + K^B^P^ < 0; 
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this inequality can be rewritten in the LMI form by recall- 
ing that a similarity transformation preserves the eigen- 
values. Therefore, defining Q = and Y = KQ, the 
closed-loop upper-bound minimization problem can be 
written a&^^ 



H 



mm 7 
AQ + QA 
I< Q < 7I, 
Q Y« 
Y m'I 



BY + Y^B^<0, Q = Q^>0 



> 0, 



where the last, additional inequality ensures a limit in 
the control effort in the form maxt>o ||u|p < n^. The 
problem pH)) has to be solved for Q, Y and 7; controller 
gains are obtained from K = YQ~^. This problem is 
again a LMI generalized eigenvalue problem, that can be 
solved using standard methods--. 



V. RESULTS AND DISCUSSION 

LMI controllers are designed wavenumber-wise, using 
the model of the system ([5]) and the design equations 
(jlOp . In particular, we consider each actuation compo- 
nent (u, f or It; on both walls) independently; a limit on 
the control effort /j, = 10, kept constant in wavenumber 
space, is used in the design of all controllers. In the de- 
sign procedure, the linear equations pertaining to each 
wavenumber pair are discretized using TV = 100 Cheby- 
shev polynomials, and modal truncation at Nt = 54 (the 
maximum affordable size of the computational problem) 
is performed prior to the actual solution of the general- 
ized eigenvalue problem (|10|) . In addition to removing 
poorly resolved dynamics, modal truncation proved to 
be necessary due to the exacting memory requirements 
of the existing LMI solvers (scaling as « N^); in perform- 
ing modal truncation, it was thoroughly verified that the 
reduced order model preserves the linear transient energy 
growth of the unactuated system. 



A. Linear analysis 

The performance of LMI controllers is evaluated first 
in the linear setting at Re = 2000; in particular, pertur- 
bations at two representative wavenumber pairs are con- 
sidered, namely, an oblique wave (a = l,/3 = 1) and a 
streamwise vortex (a — 0, P — 2). In the results reported 
here, the effectiveness of different actuation components 
is also addressed. 

The analysis reported in sec. Illll shows that the closed- 
loop system will have a non- normal behavior; it is there- 
fore natural to contrast the maximum closed-loop tran- 
sient growth G with the open-loop one, in order to verify 
that a consistent reduction in G is obtained via the min- 
imization in (1101) . Results, obtained for the same trun- 





Obliqu 
Open-loop 


e waves 
Closed-loop 


Streamwi 
Open-loop 


se vortices 
Closed-loop 


u - actuation 


67.62 


33.04 


785.93 


635.83 


V - actuation 


66.81 


11.84 


33126.72 


129.23 


w - actuation 


67.13 


32.70 


163123.13 


85.90 



TABLE I: Maximum transient energy growth G for the open- 
loop and closed-loop cases, for the oblique wave case {a = 
1,P = 1) and the streamwise vortex case (a = 0,P = 2), using 
different actuation components on both walls. Re — 2000. 



cated system used in the design and for the two wavenum- 
ber pairs considered, are reported in table HI It is shown 
that, in both the oblique wave case and the streamwise 
vortex case, solution of leads to a closed-loop sys- 
tem experiencing a reduced maximum transient energy 
growth. In particular, for the oblique wave case, v is 
the most effective actuation component, whereas actu- 
ating with w is most effective in the streamwise vortex 
case. The performance of the present controllers is fur- 
ther compared against full-state controllers designed with 
the LQR approach^'^, considered in a similar transition 
problem by Hogbcrg ct al.^. The aim of the LQR control 
is the minimization of the time integral of the perturba- 
tion energy, while keeping the time integral of the control 
effort as low as possible; in fact, the control objective is 
given in terms of the closed-loop minimization of a func- 
tional in the form 



r + oo 

Jo 



dt. 



This is substantially different from the control objective 
of the present LMI formulation pUj) . which considers 
bounds on the disturbance energy and control expendi- 
ture; therefore - at a fixed control expenditure - LMI 
controllers can be used to estimate a possible best per- 
formance (in terms of peaking suppression) of other con- 
trol strategies. In order to present a fair comparison be- 
tween the LQR and the LMI formulation, we iteratively 
design and test a LQR controller keeping Q fixed (the 
same used in LMI design) and p as a free parameter, and 
we evaluate the integral of the control energy u^u dt 
in closed loop until it matches the value computed for 
the LMI controller. The closed-loop systems, controlled 
via both LQR and LMI gains, are tested against the re- 
spective optimal perturbations, using u-actuation for the 
oblique wave case and w-actuation for the streamwise 
vortex case (the best performance cases reported in ta- 
ble |T|; further, for the streamwise vortex we consider the 
antisymmetric (with respect to the y = plane) optimal 
perturbation. The time evolution of the perturbation en- 
ergy is displayed in fig. [3] and HI Results show that, at a 
given control expenditure, the worst-case initial condition 
for the LMI-controlled system experiences a lower am- 
plification than the corresponding perturbation for the 
LQR-controlled system; the peak for the LMI-controlled 
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system occurs at later times for the streamwise vortex 
case. Despite the mild reduction in maximum amplifi- 
cation, the results shown here suggest that a control de- 
sign technique directly targeting the growth mechanism 
is able to better exploit the degrees of freedom in the 
controller to achieve a minimal transient peaking of the 
energy. 
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FIG. 3: Time evolution of the perturbation energy for the 
closed-loop system controlled using LMI gains (-) and LQR 

gains ( ), actuation using v. The initial condition for the 

two systems is the respective closed-loop optimal perturba- 
tion, and the LQR gains are iteratively designed so that the 
energy expense over the simulated time horizont for the two 
closed-loop systems is the same. Results for (a = l,/3 = 1) 
at Re = 2000. 
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FIG. 4: Time evolution of the perturbation energy for the 
closed-loop system controlled using LMI gains (-) and LQR 

gains ( ), actuation using w. The initial condition for the 

two systems is the respective closed-loop antisymmetric (with 
respect to the y = G plane) optimal perturbation, and the 
LQR gains are iteratively designed so that the energy expense 
over the simulated time horizont for the two closed-loop sys- 
tems is the same. Results for (a = 0, /? = 2) at Re = 2000. 

The performance of LMI controllers is also evaluated 



against optimal initial conditions for the unactuated flow, 
and results are reported in fig. [5] and [HI Since the 
closed-loop system has additional state equations associ- 
ated with the dynamics of wall velocity components used 
as actuators, these values are set to zero, assuming that 
at initial time the open-loop flow satisfies the no-slip and 
no-transpiration condition at the walls. In particular, for 
the oblique wave case, it is shown that actuating with v 
reduces the maximum amplification of the optimal dis- 
turbance by a factor « 8.2, whereas a less efi^ective reduc- 
tion (by a factor « 2.5) is obtained using u or w. Further, 
the growth curves in these latter cases are very close to 
each other, as the effect of actuators on the oblique wave 
is symmetric. In the antisymmetric streamwise vortex 
case, the most effective components are v and w (reduc- 
tion by a factor w 6.3 and ~ 9.9, respectively), whereas u 
has a quite poor performance (amplification reduced by a 
factor « 1.2). The differences in performance between u 
and w may be interpreted with a geometric argument, as 
for a streamwise-invariant perturbation the u component 
acts in a weakly controllable direction. It is also note- 
worthy that, in all these cases, the control action reduces 
the time interval after which the initial disturbance gets 
to its maximum amplification. 
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FIG. 5: Linear dynamics of the perturbation energy for the 
unactuated case (-), and closed- loop case actuating with u 

(-•), « ( ), w {-+)■ Results for (a = 1,^9 = 1) at 

Re = 2000; in all cases, the initial condition is the optimal 
disturbance for the unactuated flow. 

It is finally worth emphasizing that, when the closed- 
loop system experiences initial conditions in the form of 
optimal perturbations for the unactuated flow, the LMI 
and LQR performance is practically equivalent, for a 
given global control effort. Considering the best perform- 
ing LMI controllers {v and w actuation, respectively), 
results are given in fig. [7] and [51 for the oblique wave 
and streamwise vortex. For the optimal oblique wave 
(fig[7]), the linear evolution of the perturbation energy us- 
ing LQR control matches almost perfectly that obtained 
with the LMI controller. In the case of antisymmetric 



7 



500 




100 200 300 400 500 600 
t 



FIG. 6: Linear dynamics of the perturbation energy for the 
unactuated case (-), and closed- loop case actuating with u 

(-•), V ( ), w (-+). Results for (a = 0, /3 = 2) at Re = 

2000; in all cases, the initial condition is the antisymmetric 
optimal disturbance for the unactuated flow. 



streamwise vortex, we obtain a slightly larger maximum 
amplification for the LQR, that is however followed by a 
faster transient to zero if compared to the LMI case. The 
results reported in fig. [7] and E] are substantially indepen- 
dent on further decrease of the value of p: no significant 
changes in the time evolution of the perturbation energy 
are obtained, but at a far larger expense. This indicates 
that these results are close to the limit of small control 
weight. 
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FIG. 7: Comparison of the linear evolution of the pertur- 
bation energy in unactuated case (-), closed-loop using LMI 

controller (~ ■), and closed-loop using LQR controller ( ). 

Resuhs for (a = = 1) at Re = 2000 (p = 0.5), actuation 
with V, optimal perturbation for the unactuated flow used as 
initial condition. 
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FIG. 8: Comparison of the linear evolution of the pertur- 
bation energy in unactuated case (-), closed- loop using LMI 

controller (— •), and closed-loop using LQR controller ( ). 

Results for (a = 0, /3 = 2) at Re = 2000 (p = 8.0), actuation 
with w, antisymmetric optimal perturbation for the unactu- 
ated flow used as initial condition. 

B. Closed-loop transition thresholds 

After these numerical experiments in the linear set- 
ting, the performance of LMI-based controllers in terms 
of transition delay capabilities has been verified us- 
ing Direct Numerical Simulations (DNS) of transitional 
Poiseuille flow at Re = 2000, using an existing computer 
code and computing system^Si. Controllers are tested 
against initial conditions in the form of: 

• a pair of oblique waves {ao — I, (3q — ±1), in a box 
of size 27r X 2 X 27r; 

• antisymmetric streamwise vortices (ao = Oj/^o = 
2), in a box of size 27r x 2 x tt. 

These initial conditions are obtained by computing the 
optimal perturbations for the unactuated flow at the cor- 
responding wavenumber pairs, and specifically for the 
streamwise vortices the optimal initial condition having 
antisymmetric (with respect to the centerplane y — 0) 
distribution of wall-normal velocity is considered, as it 
provides the optimal transition time in the nonlinear 
case^^. Random noise, in the form of a random combina- 
tion of the first 30 Stokes modes - ordered by decreasing 
real part of the corresponding eigenvalues - is added on 
the wavenumber array (0,±l,±2)ao and (0, ±1, ±2)/3o, 
and the noise energy is chosen as 1% of the total pertur- 
bation energy. The resulting optimal perturbations are 
identical to those reported in previous works*'^^. 

In order to reduce the computational problem of con- 
trol design to an affordable size, LMI controllers are de- 
signed on the same array of wavenumber pairs where 
random noise is introduced. Furthermore, the control 
effort tuning parameter is set at ^ = 10, a value which 
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is derived from preliminary tests and previous work on 
the control of the linearized dynamics of streamwise 
vortices22i. The value of the parameter n has been kept 
constant in wavenumber space; however, it should be em- 
phasized that this parameter could be a function of the 
wavenumber pair - thus providing room for optimization 
of the control performance. 

The performance of LMI controllers is quantified by 
evaluating the closed-loop transition threshold^S of a 
given initial condition, for all the wall actuation com- 
ponents. The mixed spatial discretization (fourth-order, 
compact finite differences in y direction and Fourier ex- 
pansion in X and z directions) employs 64 grid points in 
y and 16 x 64 modes in x and z. Time integration is per- 
formed via the usual semi-implicit approach, where non- 
linear terms are advanced explicitly using a low-storage 
Runge-Kutta algorithm whereas linear diffusion terms 
are advanced implicitly via a Crank-Nicholson scheme. 
Each simulation was run over a time window of 2000 
nondimensional time units, that proved sufficiently long 
to ensure that a laminar or a turbulent state was reached 
after the initial transient growth of the perturbation en- 
ergy. In order to obtain the thresholds reported in table 
nil a bisection algorithm was employed; this procedure 
requires a large number of simulations, corresponding to 
approximately 3 months of CPU time. 

The open-loop transition thresholds reported in ta- 
ble |TI] agree with previous findings?-. Results summa- 
rized in the table indicate that the LMI controller is able 
to increase the transitional energy of the initial condi- 
tions considered. In particular, a synthetic performance 
measure is indicated in the table as improvement factor 
(I.F.), corresponding to the ratio between the threshold 
energy computed in the controlled case over that corre- 
sponding to the unactuated flow. In general, for both 
the oblique waves and the streamwise vortices, actuation 
with the wall-normal velocity v outperforms actuation 
with the other components. This behavior is expected, 
as forcing with wall-parallel components affects the flow 
by means of viscous diffusion only, whereas forcing with 
V introduces an additional non-zero momentum flux at 
the boundary. The improvement factors associated with 
the oblique wave case are higher than those pertaining 
to the streamwise vortex case, when using u and v actu- 
ation. For the u-component case, this is coherent with 
previous works^., and can be interpreted physically by 
the argument that targeting oblique waves mitigates the 
subsequent development of streamwise vortices, therefore 
reducing the strength of the associated streak instability. 
The u component provides the overall worst performance 
in the streamwise vortex case, a result that can be inter- 
preted as a consequence of the particular geometrical con- 
figuration. In fact, with respect to a streamwise-invariant 
spatial structure, u actuation works in an approximately 
null direction, whereas actuation with w is better-suited, 
as shown by its improvement factor. A similar geometri- 
cal interpretation can be given for the almost equal im- 
provement factors obtained with u and w, when an pair 





Oblique 


waves 


Streamwise vortices 




Thres. 


I.F. 


Thres. 


I.F. 


Open-loop 


2.39 ■ 10"'' 




6.47 • 10~® 




u - actuation 


9.89 ■ 10"^ 


f» 4.14 


8.57-10"'' 


f» 1.32 


V - actuation 


3.04 ■ 10"^ 


12.72 


4.86 • 10"^ 


« 7.51 


w - actuation 


9.77 ■ lO"'' 


f» 4.09 


3.86 • 10"^ 


^ 5.97 



TABLE II: Open-loop and closed-loop transition thresholds, 
as measured by DNS. The values of transitional energy are 
given with with an uncertainty of ±3%. The column labelled 
I.F. indicates the improvement factor in the closed loop with 
respect to the unactuated case with an initial condition having 
the same spatial structure. 



of oblique waves is given as initial condition. 

Results reported in table |ll] for the thresholds obtained 
using V actuation are in qualitative agreement with previ- 
ous work using the LQR controllers*; however, in quanti- 
tative terms the LQR approach outperforms the present 
LMI approach. In particular, the improvement factors 
reported with LQR for oblique waves and streamwise 
vortices at Re = 2000 are I.F. = 102 and I.F. = 10, 
respectivelyS,; therefore, LQR controllers seem to perform 
substantially better than LMI controllers in presence of 
oblique waves as initial conditions. It should be empha- 
sized, however, that such comparison is not entirely ap- 
propriate. In fact, even if the same energy norm is used to 
quantify the magnitude of velocity disturbances, control 
laws are designed with different parameters constraining 
the control effort. As a matter of fact, the linear results 
reported in fig. [5] and [6] show that the LMI controller 
performs similarly to the LQR controller for a value of 
the LQR control weight p = 0.5 and p = 8.0, respectively. 
These values are different from the value of p = 0.01, used 
uniformly in wavenumber space in Hogberg et ali^. A fair 
comparison between the two approaches is not possible in 
this case; therefore, it is impossible to draw a conclusive 
statement about the effectiveness of feedback minimiza- 
tion of transient growth versus feedback minimization of 
the disturbance energy in transition delay. 

As the ultimate goal of LMI controllers is that of pre- 
venting transition to turbulence, it can be important to 
quantify the energy efficiency of these controllers in the 
nonlinear case. For instance, given a transitional initial 
condition, it is possible to compare the energy expen- 
diture of the controller to prevent transition with the 
additional energy to be introduced into the unactuated 
flow to compensate for the increase in friction, over the 
same time window (i.e. the time necessary for the tran- 
sient in the controlled flow to die out). Referring to the 
best performing cases in table [III we consider actuation 
with wall-normal velocity v and, in the two cases, initial 
conditions having energy « 3% below the corresponding 
closed-loop threshold. Using the present nondimension- 
alization, a conservative estimate of the energy required 
for the control action in the time interval [0, T] can be 
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given bj^: 

where p is the fluctuating waU pressure and Auj the 
upper and lower blowing/suction surfaces, whereas the 
additional energy required to drive the unactuated flow 
against the increased viscous drag on the same time in- 
terval is given by: 

The ratio E^/E^ reads about 1.39 • lO'^ and 1.68 • lO'^, 
for the oblique waves and streamwise vortex case, re- 
spectively. Furthermore, linear tests using v actuation 
have shown that E^ can be of the same order of magni- 
tude of the actual reduction in maximum transient en- 
ergy growth. Hence, these results indicate that, even if 
the energy expenditure due to the control action is com- 
parable to that experienced in the linear amplification 
of the optimal disturbance, it is nevertheless negligible if 
compared to the potential energy saving due to transition 
prevention. 

VI. CONCLUSIONS 

The present work has considered the design of full-state 
feedback controllers specifically targeting the transient 
energy growth mechanism in laminar channel flow. It 



has been shown that full transpiration at both walls and 
full-state knowledge are not sufficient to ensure a mono- 
tonically stable closed-loop system via a linear feedback 
law. Further, an advanced control design technique - 
based on a LMI formulation ~ has been employed to de- 
sign feedback controllers that have been tested in both 
the linearized setting and in nonlinear, transitional flows. 

Linear tests indicated that the LMI strategy allows to 
obtain a consistent reduction of the maximum open-loop 
transient growth. At a given global control expenditure 
over a time window sufficiently long for the perturbations 
to decay to zero, the LMI-controUed closed-loop system 
experiences a lower transient energy amplification than 
a LQR-controUed closed-loop flow. However, in presence 
of an optimal perturbation for the unactuated flow, the 
performance of the two control strategies is practically 
equivalent. Results obtained in the linear setting further 
indicate that, in the case of perturbations in the form of 
oblique waves and streamwise vortices, the most effective 
actuation components are v and w, respectively. 

In the nonlinear case, it has been found that these 
controllers are capable of increasing the threshold energy 
for transition when initial conditions are given to the 
flow in the form of oblique waves or streamwise vortices; 
the effectiveness of different actuation components has 
been addressed, indicating that wall blowing/suction is 
most effective in providing a higher closed-loop thresh- 
old energy. Additionally, in transitional conditions, LMI 
controllers prove to be energy-effective, as the energy re- 
quired by the control action is negligible when compared 
to the energy saving due to avoiding transition. 
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